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Abstract 



We present a general theory for the intermediate temperature (T) properties 
of Heisenberg antiferromagnets of spin-S* ions on p-leg ladders, valid for 25";? 
even or odd. Following an earlier proposal for 2Sp even (Damle and Sachdev, 
Phys. Rev. B 57, 8307), we argue that an integrable, classical, continuum 
model of a fixed-length, 3- vector applies over an intermediate temperature 
range; this range becomes very wide for moderate and large values of 2Sp. 
The coupling constants of the effective model are known exactly in terms of 
the energy gap above the ground state A (for 2Sp even), or a crossover scale Tq 
(for 2Sp odd). Analytic and numeric results for dynamic and transport prop- 
erties are obtained, including some exact results for the spin-wave damping. 
Numerous quantitative predictions for neutron scattering and NMR experi- 
ments are made. A general discussion on the nature of T > transport in 
integrable systems is also presented: an exact solution of a toy model proves 
that diffusion can exist in integrable systems, provided proper care is taken 
in approaching the thermodynamic limit. 
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I. INTRODUCTION 



One-dimensional Heisenberg antiferromagnets are strongly-interacting quantum many 
body systems for which a detailed quantitative confrontation between theory and experi- 
ment has been possible. A rather precise, and parameter-free understanding of their low 
temperature dynamic properties has emerged in a number of recent NMR experiments per- 
formed by Takigawa and collaborators |]l|,|^,|^J^ . These systems can therefore serve as useful 
springboards towards deciphering the behavior of interacting systems of greater complexity. 

The past theoretical work on the dynamic properties of these quantum antiferromagnets 
has focused mainly only the universal behavior in the asymptotic low temperature (T) regime 
T ^ p,^,|^,|^J^. In the present paper we will extend the theory to a separate range of 
intermediate temperatures. We shall argue that under suitable conditions, to be described 
precisely below, this intermediate temperature range can be quite wide, and is described by a 
continuum dynamical model quite different from that required for T — * 0. This intermediate 
temperature dynamics was discussed briefly for a limited class of antiferromagnets in the 
last section of Ref . 

Our work will also connect with earlier investigations of the dynamics of classical lattice 
antiferromagnets (Refs P,p!0| and references therein). In a sense, our paper provides a bridge 



between the modern quantum dynamics and the classical studies of the 70's. There is an 
overlapping window of validity for our theory and the classical investigation of Reiter and 
Sjolander |jlO|, and here our results are generally consistent with theirs, although there are 



some details that disagree. We will review this earlier work, in the context of our results, in 
Section 0. 

A large fraction of the experimental examples of one-dimensional Heisenberg antiferro- 
magnets consist of p parallel, coupled chains of spin S ions (for p = 1 these are ordinary spin 
chains, while for p > 1 these are commonly referred to as p-\eg ladders). For all T < T^^xj 
where T^l^ will be defined shortly, these antiferromagnets are described by a universal quan- 
tum field theory: the one- dimensional 0(3) non-linear sigma model. This field theory has 
the quantum partition function (in units with h = ks = 1, which we use throughout) 



Zq = J V'n{x, r)5(n^ — 1) exp ^— J d-^ J dr C 



2cg 



zH X n I + 



dx I 




Here n(a;, r) is a three-component unit vector representing the orientation of the antifer- 
romagnetic order parameter at spatial position x and imaginary time r, c is a spin-wave 
velocity, and H is a uniform external magnetic field — we will be interested only in the linear 
response to H. There are two dimensionless coupling constants in £, 9 and g. The first, 9, 
is the coefficient of a topological term, and has the value 6' = vr for 25*^ an odd integer, and 
the spectrum of excitations above the ground state is then gapless. For 25*^ an even integer, 
6' = 0, and then there is gap to all excitations. The coupling g plays a role in determining 
the energy scale at which certain crossovers (to be discussed below) take place, but does not 
modify the physics otherwise. A straightforward semiclassical (large S) derivation shows 
that 
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FIG. 1. Renormalization group flows for the dimensionless coupling g in (IJ-)- -^^^ ^ = 0; 5 
a runaway flow to 5 = oo and the ground state is a quantum paramagnet with a gap A. For 9 = n, 
there is a fixed point at g = gc, of order unity, and near it the flow is dg/di oc {g — gc)"^. This fixed 
point is described by the /c = 1, SU{2) Wess-Zumino-Witten model. The crossover between the 
g = and g = gc fixed points takes place at an energy scale of order Tq. The region g > gc usually 
corresponds to a gapped state with spin-Peierls order, and is not considered in this paper. 
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where J is the exchange constant along the legs of the ladder, J± is exchange on the rungs, 
and a is the lattice spacing along the legs; we have assumed here a model with only nearest 
neighbor exchange, but the estimates g ~ 1/Sp and c ~ JSa hold far more generally. An 
important observation is that g becomes small for either large S or p. We will be especially 
interested in the small g case in this paper. 

Let us now discuss the value of T^^^x below which ( |1.1|) holds. The basic argument follows 
that made by Elstner et al [|lT| m d = 2. At a temperature T, the characteristic excited 
spin- wave has wavelength c/T, and the continuum quantum theory will apply as long as 
this wavelength is longer than the lattice spacing, a, of the underlying antiferromagnet. For 
p-leg ladders, description by a one- dimensional quantum model requires that the wavelength 
be larger than the width of the ladder, pa [0]. Using the value of c in ( |1.2| ), our estimate 
for T^}^ is then 



2JS 
p 




.1/2 
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To reiterate, the quantum theory ( |1 . 1| ) applies to the lattice antiferromagnet at all T below 
that in ( 101) . 

Let us now review the well-known, T = 0, renormalization group properties of ( |1 . 1| ) ||13|| . 
The topological angle 9 remains fixed at 6' = 0, vr, while the flows of the coupling g are 
sketched in Fig |I|. For both cases ^ = 0, vr there is a fixed point (7 = which is unstable at 
low energies. Indeed the beta function describing the flow away from (7 = is independent 
of 6 to all orders in g. However, non-perturbative, topological effects do distinguish the two 
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values of 6. For ^ = 0, the flow is believed to continue all the way to g = oo, corresponding 
to a quantum paramagnetic ground state with an energy gap A. In contrast, for 6 = n, the 
flow is into a strong coupling, infrared stable fixed point ai g = gc- There is a scale-invariant 
and gap less theory which describes this fixed point — the k = 1, SU{2), Wess-Zumino-Witten 
model. For g > gc, there is a runaway flow to g = oo, usually associated with the appearance 
of spin-Peierls order; this last regime will not be discussed in this paper. 

Our primary interest here shall be in the region in the vicinity of the unstable g = 
fixed point. For both 6* = 0, tt, there is a characteristic energy scale, usually denoted Ajjg, 
which determines the location of the crossover from the vicinity of the g = fixed point to 
the strong-coupling behavior. For energies or temperatures smaller than A;^, the strong 
coupling behavior should apply, and as we have just discussed, this is quite different for 6* = 
and 6 = 71. However, the physics at energies or temperatures larger than is controlled 
by the flow in the vicinity of the g = fixed point, and this is common to both ^ = 
and 9 = n. The can be estimated from the structure of the perturbative 

/3-function. For small g 



A^ ~ Jexp 1^-- -J ; (1.4) 
we have neglected here a prefactor of a power of g coming from higher-loop corrections. So 



from ( |1.2|) , for either S or p moderately large, the scale Aj^s becomes exponentially small. 

We are now ready to discuss the static thermodynamic properties of the quantum field 
theory (|1.1|) as a function of T. We will characterize the system by T dependence of two 
important observables: Xu(T) and ^(T). The first is the uniform susceptibility, Xu, which is 
the linear response to the field H = {0,0, H): Xu = (T / L){(P In ZQ/dH^)\H=o, where L is 
the (infinite) length of the spatial direction; this is the susceptibility per rung of the ladder. 
The second is the correlation length, ^(T), which determines the exponential decay of the 
equal-time two-point n field correlator as a function of x. 

We will consider temperatures above and below Aj^g in turn. 
(A) T<Aj^ 

For T < Ajjg, as just noted, we must distinguish 9 = and 9 = n. 
(A.l) 9 = 

For ^ = 0, there is an energy gap A, and the susceptibility is simply that of a dilute, 
thermally activated, classical gas of triplet magnons above the gap; these contribute an 
exponentially small susceptibility ||l^ , p!5| 



/ 2A \ 

Xu{T) = [-^) e-^/^ ; T<Aj^,9 = 0. (1.5) 



Experimentally, we can view ( |1.5|) as the definition of the gap A and the velocity c, which 



are to be determined by fitting measurements to ( |1.5| ). The correlation length, ^(T) takes 
a finite, T-independent value in this quantum paramagnet, up to correlations exponentially 
small in A/T; for the case where A is significantly smaller than J, we have [P^,P^ 

= 1 ; T<Aj^,9 = 0. (1.6) 
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(A.ii) e = Ti 

For the gapless case, 6 = ir, there are excitations with non-zero spin at arbitrarily low 
energies, and so Xu{T) remains non-zero as T ^ p!8|JT9|] : 



Xu{T) 



2-KC 



1 ln(ln(ro/r)) 
21n(To/T) 41n2(To/T) 



+ . . . 



T<^Ms^e 



TX. 



(1.7) 



Again, this experimentally defines c and a new temperature scale Tq which determines the 
onset of a logarithmic correction to the T = susceptibility due to the slow fiow into the 
fixed point ai g = g^. The correlation length of this critical paramagnet now diverges as 



c 



1 + 



1 ln(ln(To/T)) - 

21n(ro/r) A\n\To/T) 



+ . 



T<Aj^,e 



TT. 



:i.8) 



The dynamical properties of quantum antiferromagnets in the low temperature regime T < 
have been discussed at length in Ref for the gapped case (6* = 0) and in Refs [§|,H,|3)|H 
for the gapless case {9 = n). 

Note that both A and Tq are energy scales characterizing the fiow into the strong-coupling 
region. These should therefore be universally related to Ap-^ which is the scale of fiows out 
of the weak coupling region. We will discuss the universal relation shortly, once we have 
defined more precisely. 
(B) T > A^ 

Let us now consider the regime Ajjg < T < Tj^^^^; the upper-bound is necessary to ensure 
the continuum quantum theory still applies. The existence of this intermediate temperature 
regime requires that < T^l^, a condition that is not well satisfied for small S and p, 
and so this regime almost certainly does not exist for p = 1 and S = 1/2, 1, but there is 

r7| . Here we are controlled by physics in the 
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evidence that it is present for p 
vicinity of the g = fixed point, and it should be possible to treat quantum fiuctuations in a 
renormalized perturbation theory in g. As discussed in Ref Q], a non-perturbative treatment 
of the thermal fiuctuations is still necessary, but this can be carried out exactly because of 
the low spatial dimensionality. The result of such a calculation is M 
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where 7 is Euler's constant. This result can also be viewed as the precise experimental 
definition of A;^. For completeness, we also quote the result in this regime for the correlation 
length, e(T): 



e(T) 



2nT 



47re~^r T 
ln( — : ) + Inln 



A 



MS 



A 
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MS 



^Ais < T < T^l^ 
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We now give the promised universal relationship between the energy scales characteriz- 
ing the weak {Ajj^) and strong (A, Tq) coupling regimes. This requires non-perturbative 
knowledge of the renormalization group fiows, and can only be obtained from an analysis of 
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the full thermodynamic Bethe-ansatz solution of the quantum field theory For ^ = 0, 

such an analysis was carried out in Refs [21|, and the result is now well known: 



A- 



MS 



6 = 0. 



(1.11) 



For the gapless case, 6 = n, we will present a derivation of the required relationship in 
Appendix A, building upon some recent results p3| , |22[| ; our result is: 
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MS 
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TT. 



(1.12) 



We are now finally in a position to state precisely the main objective of this paper. We will 
describe the dynamical properties of one-dimensional antiferromagnets in the intermediate 
temperature regime Aj^g < T < Tj^l^. We quickly note as an aside that most of our results 
actually hold over a wider regime of temperatures Ajjg <T < T^l^, where we will define the 
discuss the origin of T^l^ > T^^^x below; for now we ignore this point. A common treatment is 
possible for the 6 = and 6' = vr in this regime, with the two cases differing only in the input 
values of the static parameters Xu{T) and ^(T) as defined by (|1.9| , |1.10| , |1.11| , |1.12D . Moreover, 



as Ajj^ becomes exponentially small for moderate values of S or p, this regime can be quite 
wide and should be readily observable experimentally. Indeed, there is good evidence from 
recent measurements of static properties in quantum Monte Carlo simulations |]T^ that this 
universal intermediate temperature regime exists even for 5* = 2 spin chains. 

The formulation of the dynamics properties for A^g^ < T < Tj-^^^ was already discussed 
in Ref. [Q. The key point ||^ is to notice that the energy of a typical spin- wave excitation, 
which is of order c^~^{T), is parametrically smaller than T when ^ obeys ( |1.10D . So the 
thermal occupation number of these spin-wave modes is large: 



T 



> 1. 



:i.i3) 



The second expression in (|1.13|) is the classical equipartition value, which indicates that the 
spin-wave excitations may be treated classically. The classical partition function controlling 
these fluctuations can be deduced by demanding that its correlations match with (l.£, 1.10 ), 
while the dynamic equation of motions following by replacing the quantum commutators 
with Poisson brackets. In this manner, the problem reduces to the effective classical phase- 
space partition function M 



Zc = j Vn{x)Vl.{x)5{n^ - 1)5(L ■ n) exp 



c 



T 



n 



c 



dx 



T^(T) 



dx 



+ 



Xu±{T)' 



;i.l4) 



Here n(x) is a classical variable representing the orientation of the antiferromagnetic order 
and L(x) is its classical, canonically conjugate angular momentum. Because n(x) is of unit 
length, its motion is always in a direction orthogonal to its instantaneous direction, and 
there is no radial kinetic energy: the square of angular momentum represents the entire 
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kinetic energy, and Xu± is the moment of inertia of the fluctuating n. The value of Xu,± 
can be determined by realizing 0] that L(x) is simply the classical limit of the quantum 
operator corresponding to the magnetization density 

6C 

L=-^; (1.15) 

then demanding that Zc reproduce the correct uniform susceptibility to a uniform external 
field under which He — ^ He — J dxH ■ L, we obtain 

XuiT) = ^Xu±{T)- (1.16) 

the factor 2/3 comes from the constraint L ■ n = 0, so there are only two independent 
components of L at each spatial point. 

Notice that ( |1.14| ) involves a functional integral over the commuting fields n and its 



conjugate momentum L as functions only of the spatial co-ordinate x, but is independent 
of the real time, t. It therefore yields only equal-time correlation functions, as is the usual 
situation in classical statistical mechanics. To obtain unequal time correlators, we have to 
separately specify the equations of motion, and these are obtained by replacing quantum 
commutators with Poisson brackets. For the fields n{x), L(x) these are 

{Laix), Lf3{x')}p^ = e«/3^L^(x)(5(x - x') 
{La{x),nf^{x')}pj^ = eap^n^{x)6{x - x') 

{na{x),np{x')}p^ = 0. (1.17) 

where a, (3'y = 1, 2, 3. The equations of motion (in real time) now follow from the Hamilto- 
nian He, and they are 

On 1 ^ 

dt Xu±{T) 

|^™„x0 (1,18) 



We are now interested in unequal time correlation functions of ( 1.18| ), averaged over the 
classical ensemble of initial conditions specified by Zc- 

To complete the quantum-to-classical mapping, we recall 0] the relationship between 
the correlations of the underlying antiferromagnet and the quantum field theory Zq, and 
those of the classical non- linear wave problem defined by ( [1.14| ) and (|1.18|) . Correlations 



of the antiferromagnet in the vicinity of the antiferromagnetic wavevector are given by the 
correlations of n(x,t) under Zq., and are related to those in the classical problem by @] 



(n(x,t)-n(0,0))Q = ^ 



(n(x,t)-n(0,0))c, (1.19) 



where the subscript Q represents averages under the quantum partition function Zq, the 
subscript C represents averages under the classical dynamical problem defined by ( |1.14| , |1.18D , 



7 



and A is an overall T-independent normalization related to the amplitude of the correlations 
at T = 0. Next, correlations of the antiferromagnet in the vicinity of zero wavevector are 
given by the correlations of L(x,t) under Zq, and are essentially equal to those in the 
classical problem by 0] 

(L(x, t) ■ L(0, 0))q = (L(x, t) ■ L(0, 0))c. (1.20) 

The absence of an overall rescaling factor here is related to the conservation of the total 
magnetization density. 

The main objective of this paper is to evaluate {n{x,t)n{0,0))c and (L(a;, t)L(0, 0))c- 
An important property of these correlators is that they satisfy simple scaling laws which 
allow us to completely scale away all dependencies on ^(T) and Xu±iT), and to express 
everything in terms of parameter-free, universal functions. These scaling laws follow from 
the fact that ( |1. 14 , 1.18 ) define a continuum classical problem which is free of all ultraviolet 



divergences: this will become evident from our analytic computations in Section |I| and 
the numerical results of Section |T|. Consequently, simple engineering dimensional analysis 
involving rescaling of x, t and L can be used to absorb dependencies on the dimensionful 
parameters. In this manner, it is not difficult to show that M 



(n(x,t)-n(0,0))c = $„ (x,t) 

^(T) 



(L(x, t) ■ L(0, 0))c=( <^>L {x, t) , (1.21) 



where 



_ X 
X = 



t = t 



T 



.aT)Xu±iT) 



1/2 

(1.22) 



and $n(^, t) and $L(a;, t) are universal scaling functions. It is the primary task of this paper 
to determine these scaling functions. In principle, these scaling functions are determined 
by solving ( |1.14|jl.l8| ) after setting all parameters equal to unity, T = ^ = Xu± = 1, while 
replacing x,t by x,t. Notice that the resulting equation is then parameter- free, and so there is 
no explicit small parameter in which any kind of expansion can be carried out; nothing short 
of an exact solution will do. We think this reasoning invalidates some of the conjectures 
on exactness of results made in Ref [TU|, as we will discuss further in Section |V|. However, 



when the arguments of the scaling functions are themselves small, i.e. \x\, \t\ -C 1, then a 
systematic perturbation expansion is possible, and this will be presented in Section |T|. 

We close this introductory discussion by returning to the issue of the maximum temper- 
ature up to which these results can be applied to lattice antiferromagnets. The appearance 
here of a classical spin model suggests that one should think about classical spin models 
obtained by starting directly from the lattice quantum spin model, without the use of the 



quantum field theory Zq as an intermediate step [|ri|]. Such a classical description will 
only work for large S, and we can ask the question of when the resulting classical model 
can be described by a continuum classical theory, which will clearly be the one defined by 
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( |1.14| , p!rT8D . The correlation length of a classical spin antiferromagnet is of order JS'^pa/T, 
and a continuum description will work provided this is larger than the lattice spacing a and 



the width of a ladder system pa. This gives us the estimate |TT 



:i-23) 



In the regime T^l^ < T < T^l^, we can use a purely classical description of the spin 
model: its correlation length and uniform susceptibility will not be universal, but has to 
computed for the specific model under consideration. For the model with only nearest 
neighbor exchange, a standard computation on the classical antiferromagnet gives for the 
uniform susceptibility 



Xu{T) 



P 
6Ja 



1+1 




7^(1) ^ T < T*-^) 

max ■ ■ max ' 



and for the correlation length. 



JS^pa 
T 



■ T^^'^ < T < T^^^ 

5 max ^ ^ max ' 



:i.24) 



:i.25) 



It is satisfying to note that there is a precise agreement between the results ( p^|1.10|) and 
( |1.24| , pr25|) at the common boundary of their respective regions of validity, T ~ ^^ax- 
the nearest neighbor model under consideration we use the estimates for g and c in ( p..2D , 
and then using ln(T/Ajg^) ^ 27r/(yf, we find the required agreement. 

It appears useful to review the final status of the regimes of validity of the model studied 
here. The universal, continuum, classical model ( 1.14 , LT^ ) describes all one-dimensional 
Heisenberg antiferromagnets in the temperature regime A^g^ < T < T^l^] this regime is 
wide and well defined for moderately large values of 5* or p, and there is evidence that it 
exists already for p = 1, 5 = 2 []T^. The complete definition of this classical dynamical 
model requires the input of the temperature-dependent static parameters ^(T) and Xu{T)- 
In the regime < T < T^l^ these parameters are universally specified by (|1.9|J1.10| )i with 
^Ws gi'^sn by for gapped spin chains {9 = 0) and by ( |1.12| ) for gapless spin chains 

{9 = 7i). In the higher temperature regime T^l^ < T < T^^x? these parameters are given 
by ( |1.24| ) and ( |1.25| ) for the model with only nearest neighbor exchange, and by related 
non-universal expressions for other antiferromagnets. 

The following sections contain technical details towards the determination of the scaling 
functions $t 



and $L in (1.21 



along with some theoretical analysis on the relationship 
between integrability and diffusion. We will begin in Section || by describing the analytical 
short time expansion of the correlators $„ and The long time limit will then be studied 



numerically in Section |ITI|. The subsequent Section |^ discusses issues which are somewhat 
peripheral to the main focus of this paper: the continuum equations of motion ( |1.18| ) are 
known to be integrable p4i, and this raises numerous fundamental questions on the nature 



of spin transport in integrable systems. These will be addressed in Section ^ by the study 
of a simple, integrable, toy model whose spin correlators can be determined in close form. 
Further, we will see that these correlators have a striking similarity to those of ( |1.18D . 
Readers interested primarily in spin chains can omit Section |^ and skip ahead to Section 
where we will describe the implications of our results for experiments. 
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II. SHORT TIME EXPANSION 



This section will determine the small t expansion of the scaling functions $n and $l 
in ( p. .211 ). Normally, there is a completely straightforward way of determining the short 
time expansion of interacting systems pSf — it can be related, order-by-order, to equal-time 
correlators involving higher moments of the fields. However, this standard procedure does 
not work for the model (|1.14 , LT8|) of interest here. This is because we are dealing with 
a continuum model with an infinite number of degrees of freedom, n(x), L(x), present at 
arbitrary short distance scales. If we naively generate the moment expansion, we find that 
the terms quickly acquire rather severe ultraviolet divergences. 

A separate theoretical tool is necessary to generate the short time expansion, and this 
shall be described here. We shall use an analog of the field-theoretic method known as chiral 
perturbation theory. As we shall see below, the expansion is actually in powers of \t\ — this 
implies a non-analyticity in the t-dependence at t = 0, which is in fact the reason for the 
ultraviolet divergences in the moment expansion. The latter method only gives an analytic 
expansion in t, by construction. 

In this section, and in Appendix B, we will use units in which T = Xu±{T) = 1. However, 
we will retain explicit dependence on ^ = ^(T). It turns out to be quite useful to keep track 
of powers of ^. Indeed, our computations will be designed to generate an expansion of the 
correlators in powers of 1/^, and this is a posteriori seen to be a short time expansion. We 
will return to physical units in stating our final results. 

We will therefore consider the problem of unequal time correlation functions of the non- 
linear partial differential equations 



dn 

dt dx"^ 



(2.1) 



when averaged over the ensemble of initial conditions defined by the partition function 



2c = J Vn{x)VL{x)6{n'^ - 1)5(L ■ n) exp ^-^ J 



dx 



dx 



+ - 2^m^n^ 



(2.2) 



The last term in the action represents a field of strength turned on in the z direction 
and serves as a regularization parameter for our perturbation expansion. At the end of our 
calculations we shall let m — > 0. From this partition function, one can immediately find the 
equal time correlation functions to be 



(L(a;,0)-L(0,0)) 
(n(x,0) -11(0,0)) 



26{x) 



(2.3) 



the subscript C is implied on all averages in this section, unless stated otherwise. 

We first construct our perturbation expansion in powers of 1/^ for the equal time problem 
specified by ( p.2|) , and check that we do arrive at the correct correlation functions as specified 
in (p. 31). The extension to the unequal time problem will then be straightforward. First, the 
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constraints on the fields n and L are solved by introducing two complex scalar fields and 



1 



Uy 

Lfr. 



2C 
1 



1 - 2(j)(f)*/^ 



iV2 



(2.4) 



We introduce this decomposition to the functional integral (|2.2| ) and expand the square 
roots in power series of 1/^. In this manner, we arrive at an interacting field theory with an 
infinite number of interactions, with 1/.^ as the small coupling. However, to any particular 
order in perturbation theory in 1/^, we only need to keep a finite set of interactions terms. We 
evaluate the correlation functions in real space using the ordinary machinery of diagrammatic 
perturbation theory. It is well known that such a perturbation expansion is plagued by 
infrared divergences. The decomposition ( |2.4| ) takes it for granted that the n field is ordered 
in the z direction and represents small spin wave fluctuations around the ordered state. 
The infrared divergences are a signature of the fact that this assumption is wrong in one 
dimension. By introducing the external field ^m^, we introduce long range order into the 
system and thus regularize the divergences. The divergences show up as poles in 1/m. But 
if we calculate the 0(3) invariant correlation functions as in (|2.3|) , we find that all poles 
cancel. Keeping this in mind, it is not difficult to show that the results of perturbation 
theory agree with ( p.3| ) at every order. The details are to be found in Appendix B. 

Once we have the correlation functions in the equal time ensemble, we proceed to evaluate 
the unequal time correlation functions. To do that we again insert ( |2.4|) into the equations 
of motion ( ^TID and expand in powers of 1/,^ to get 



'dt 



1/2 



d'^d) 1 d 
— - H 




dx J 



+ . . . 



(2.5) 



Higher order terms add further nonlinear interactions. We solve the initial value problem 
for each of the fields by an iterative strategy. First the free wave equation is solved and 
the solution is plugged back into the lowest order nonlinear term to solve the problem to 
the first order. To evaluate the correlation functions, we just multiply the fields and carry 
out the average over initial conditions. The initial condition averages are, of course, known 
from the calculations described above. The technical details are relegated to Appendix B 
and here we only quote the results to one loop order: 
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Y^jfyj (L(^, t) ■ L(0, 0)) = -t) + 5{x + ^)) (l " y) 

+^ {e{x + t)- e{x - 1)) + 



(2.6) 



(n(x,t) -11(0,0)) 



c = 1 - -(|3^+t| + \x-t\) 

+— f3(x + t)^ + 3(x-t)^ + 2|x + t||x-t|) +0(x,t)^ (2.7) 
16 ^ ^ 



where x and t are defined in (|1.22|) , and the deha function in ( p.6|) is interpreted as 0(x, t)~^ 
In these resuhs if we set t = 0, we immediately recover the equal time results ( p.3| ) to 
the corresponding order in 1/^. The structure of the correlation function (|2.6|) refiects the 
causal propagation of the conserved angular momentum L. The first term simply represents 
the free propagation of the angular momentum density which is completely concentrated on 
the "light cone". However, interactions to order 1/^ do modify the free propagation and 
transfers the angular momentum density from the surface of the light cone to its interior. 
This is represented by the second term which is non- vanishing only inside the light cone. It 



can be checked that the spatial integral of ( |2.6| ) remains independent of t, as must be the 
case due to conservation of total angular momentum. 

The result ( |2.7| ) has significant implications for the dynamic structure factor, Sik^uS)^ of 
the antiferromagnetic order parameter: 



S{k, uS) 



dxdt{n{x, t) ■ n(0, 0))Qe 



-i(kx—LL!t) 



(2.8) 



First, we note that the result (|1.19|) and the scaling form (|1.21|) , imply that the dynamic 
structure factor satisfies 



T 



1/2 



Sik)^s{k,cJ) 



(2.9) 



where S{k) is the equal time structure factor, which is known exactly (apart from the overall 
normalization A) 



S{k) 



27r 



A 



In 



S{k, u) 



T 



l + k^^^{T) 



(2.10) 



and $5 is a universal scaling function of 



k = ki{T) 



UJ = UJ 



(2.11) 
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which describes the relaxation of the equal time correlations. The prefactor in ( p.9|) has 
been chosen so that the frequency integral over $5 is normalized to unity for every k 

0° dZJ , , , 

—^sik,iJ) = l. (2.12) 

We will now show that, modulo some very mild assumptions, the result ( |2.7| ) exactly fixes 
the form of ^s{k,uj) for \k\ ^1. In this short distance regime, we are at distances shorter 
than the correlation length, and the system should look almost ordered. So we may expect 
that the spectrum consists of weakly damped spin waves, and this motivates the following 
ansatz for $5 in the regime |fc| 3> 1 and |a;| ~ \k\: 

where 7(fc) is the unknown spin-wave damping parameter; we will shortly determine the 
large k limit of 7(A;). We now have to take the Fourier transform of ( p.9[ ), ( p.lO| ) and ( p.l3[ ). 



compare the result with (|2.7|) . In making this comparison, we should keep in mind that the 
large k,uj behavior of S{k,u)) can only determine the non-analytic terms in the small x,t 
expansion. First, the integral over frequencies can be performed exactly for the form ( p.l3| ), 
and ( |2.9| ) implies 

(n(..«).n(0.0)>,^/g ''^-:!'"-'^' +... (2.14) 

k +1 

2 

We reiterate that (|2.13| ) is valid only for ^ 1, and so only the 1/k contribution from the 
{k + 1) denominator in ( p.l4| ) can be taken seriously. Also, we already have a non-analytic 
\t\ dependence in the exponential, and so we can expand this in powers of in this manner 
we reduce ( p.l4| ) to 



, , , 2 ^-rcos(kx) cos(kt) 2\t\ f°° ^-r'f ( k) cos( kx) cos( kt) 

(n(x,t) ■n(0,0))c = - / dk — _ ^ / dk^^ ^ + ..., 

vr Jfci 71 Jk2 k 

(2.15) 

where ki^2 are some large positive constants, and we have assumed that 7 is an even function 
of k. The values of the integrals over k surely depend upon p and p\ but the key observation 
is that the non-analytic terms in x and t do not. This follows from the result 

-— cos(A;a;) = — , ' ' — + ..., 2.16) 

k^ ^ ' 2r(z/) cos(7rz//2) ' ^ ' 

where all omitted terms can be written as a series in non-negative, even integer powers of x 
(this allows an additive constant, independent of x). We now assume 7(/c 00) ~ A;", and 
then demand consistency between ( p.l5| , pTT6D and the non-analytic terms in (p.7|) . It is easy 
to see that we must have a = 0, and so 7(00) is a constant. Further, the unknown additive 
constant associated with the second integral in (|2.15|) must be such that there is no single |t| 
term in the correlator. Applying (|2.16|) to (|2.15| ) with this understanding, we deduce that 
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(n(x,t)-n(0,0))c 



-(|x + t| + |x-t|)(l-7(oo)|t|) + ..., 



(2.17) 



where all omitted terms are either analytic in x and t, or involve subleading non-analyticities. 
We should now compare ( |2.17D with ( |2.7| ) by matching only the non-analytic terms in the 
vicinity of the light cone x = iit. In this latter region we can approximate \x + t\\x — t\ in 
( ^.71 ) by 2|t||x±t| — then the non-analytic terms in ( |2.17D and ( |2.7D match perfectly, and we 
obtain one of our important exact results 



7(00 J 



1 
2' 



(2.18) 



Returning to physical units via ( p. Ill) , we conclude that the frequency r(T) = (u;/ti;)7(oo) 
given by 



r(r) = - 



nl/2 



.aT)Xu±iT) 



(2.19) 



describes damping of spin- waves for \k\^(T) ^ 1. We will discuss the experimental implica- 
tions of this result in Section |V|. 

Results for the spin-wave damping have been obtained earlier by Reiter and 
Sjolander [jlO[, in their studies of classical lattice antiferromagnets. If we insert the clas- 
sical values ( |1.24| , |1.25D into (|2.19|) , we obtain 



r(T) 



T 

Sp 









1 + 


(-3 











1/2 



(2.20) 



which agrees with their result for p = 1. Keep in mind, though, that our result ( 2.19 ) 
has a much wider regime of applicability, beyond temperatures in which a purely classical 
thermodynamics holds. Further discussion on the relationship between our and earlier results 
appears in Section |V[ 



III. NUMERICAL RESULTS 

The previous section allowed us to determine the two-point correlators for small \t\. 
Here we will present numerical simulations which examine the large \t\ limit. These were 
performed on a discrete lattice realization of (|1.14|jl.l^), with lattice spacings of ^/16 or 



larger, and had only nearest neighbor couplings between the n vectors. Initial states were 
generated by thermalizing the system by the Wolff algorithm [^. The time evolution was 
carried out by a fourth order predictor-corrector method, and its accuracy was tested by 
keeping track of the conserved total energy and the lengths of the n vectors. 

Our simulations are similar to many earlier studies of classical spin chains (see Refs P7|PB 



and references therein). However, there is an important difference in that we are dealing 
with rotor variables n, L, rather than classical spins S which obey Poisson bracket relations 
like those for L in (|1.17|) . 

We will consider correlators of n and L in the following two subsections. 
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FIG. 2. The correlator (n(0,t) • n(0,0))(7 as a function oft, which is defined in ( 1.22| ). Results 
with lattice spacing ^/16 and (,/8 are shown, and their near perfect overlap indicates we have 
reached the continuum limit. 



A. Correlations of n 

We first consider dynamic correlators of n. The aim of our simulations is to obtain results 
for the dynamic structure factor, S{k, u), in regimes beyond the case \k\ 3> 1, \uj\ ~ \k\ which 
was studied by the short time expansion. Our results were obtained for two cases — at equal 
positions (local) and at zero wavevector. 

The local correlator is measured in NMR experiments, and we computed the local dy- 
namic structure factor, Si{uj) defined by 



rft(n(0,t) ■n(0,0))Qe'"* 

J-oo ZTT 

By ( |1.19| ) and ( p,.21| ), 5*^ satisfies the scaling form 



(3.1) 



S,(cu) = A 





2 








T 



1/2 



(3.2) 



where is a fully universal function (with no arbitrariness in its overall amplitude or the 
scale of its argument) with a unit integral over frequency 



(3.3) 



Our results for the local time-dependent correlations are shown in Fig ^ and its Fourier 
transform to frequency in Fig |. Two different lattice spacings were used, and the good 
overlap of the data confirms that we are examining the continuum limit. The correlations 
decay rapidly in time, but also show a brief, but clear oscillation; this oscillation results in 
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FIG. 3. The Fourier transform of Fig into frequency. This yields the universal scaling function 



$z defined in (13^ 



a finite frequency peak in Si{uj). We will discuss the physical origin of this oscillation after 
we have considered the zero momentum correlator. 

Turning to the zero momentum correlator, we express our results in the scaling form 
( ^IqD and obtain values for the scaling function ^^(OjCl;). We emphasize that the ansatz 
( p.l3| ) does not hold for k = 0. Our results for the k = correlators of n are shown in 
the time domain in Fig | and after the Fourier transform to frequency in Fig |^. As with 
the X = correlations above, the k = correlator shows a rapid decay, along with a brief 
oscillation; the latter leads to a finite frequency shoulder in $5(0, cJ). 

How do we understand this finite oscillation frequency observed in both the x = and 
k = correlators of n ? One way is to compare with the exactly known results p9| , p0| of the 
model with an A^-component vector n in the limit of large N. At N = 00, S{0, u) consists of 
a delta function at a finite frequency 00 ~ T/ ln(T/A;^). So we can view the finite frequency 
as a remnant of the N = 00 response at A^ = 3. However, there is a related, more physical, 
way to interpret it. The underlying degree of freedom have a fixed amplitude, with |n| = 1. 
However, correlations of n decay exponentially on a length scale ^(T) — so if we imagine 
coarse-graining out to ^(T), it is reasonable to expect significant amplitude fluctuations in 
the coarse-grained field, which we call 0q,. On a length scale of order ^(T), we expect the 
effective potential controlling fluctuations of (pa to have minimum at a non-zero value of 
|0Q,|, but to also allow fluctuations in about this minimum. The finite frequency in 
Figs 1^, 1^ is due to the harmonic oscillations of 0q, about this potential minimum, while the 
dominant peak at a; = is due to angular fluctuations along the zero energy contour in the 
effective potential. This is interpretation is also consistent with the large A^ limit, in which 
we freely integrate over all components of n, and so angular and amplitude fluctuations are 
not distinguished. 
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FIG. 4. The correlator / dx{n{x, t) • n(0, 0))c/C{'^) & function of t, which is defined in ( |1.22 ). 
We used a lattice spacing ^/8. By (|2.1C| ), the t = value of this should be 2, and the difference is 
due to the finite lattice spacing. 




1 _ 2 

CO 



FIG. 5. The Fourier transform of Fig ^ into frequency, which when combined with ( ^.91) , leads 
to the the universal scaling function $5(0, w). The full line is at lattice spacing ^/8, and the dashed 
line is for ^/16. The data at ^/16 is 'noisier' because of insufficient averaging for the longer time 
data. 
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FIG. 6. Numerical results (full line) for the {^{T)/Txu±{T)){L{0,t) ■ L(0, 0)) correlation func- 
tion for short times t < 1 on a lattice with spacing ^/8. The results are compared with short time 
expansion (dashed line) in (|3.6D , valid for t ^ 1. 



B. Correlations of L 



We obtained numerical results only for the a; = correlator of L. The short time behavior 
of this is given in ( |2.6| ). At long times, we expect the conservation of total L to be crucial 
in determining its asymptotic form. In particular, one natural assumption is that the long 
time correlators of C are diffusive; in this case we expect 



at large t. Consistency of this with the scaling form ( |1.21| ), implies that the diffusion 
constant, D must obey 

D = B ^^7^ (3.5) 

where ;B is a dimensionless universal number. 

Our numerical analysis of the autocorrelation was carried out on a system of 800 sites. 
The predictor-corrector method turns out to exactly conserve angular momentum and we 
maintained energy conservation to 4 significant digits over the duration of the simulation. 
We averaged over 9,600 initial conditions. 

First we tested our results against the known exact short time expansion. This is shown 
in Fig p. At these short times the lattice corrections are quite significant, and our comparison 
in Fig ^ is with the chiral perturbation theory carried out in the presence of a lattice — the 
generalization of the result ( p.6|) to a lattice model with nearest neighbor couplings and 
lattice spacing e.^ is 
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FIG. 7. Numerical results for {^{T)/Txu±{T)){L{0,t) ■ L(0,0)) as a function of t in a log-log 
plot. The smoother line is for lattice spacing ,^/8, while the noisier line is for lattice spacing 
^/16. The agreement of the two results is evidence that this decay is a property of the continuum 
limit. A straight line fit (not shown) to the ^/8 data is almost perfect, and its slope indicates 
that the correlator decays as t ' . An equally good fit to the data was obtained by the function 



ai/Vt + a2/t, with the second subleading term contributing only about a 10% correction at the 
largest t. 



(L(x,t)-L(0,0)) = 2 / e'^-cos{ukt) 

'Tv/t Zn 



1 + 



dp cos{uLipt)e^P^ 



+ 2 



r^' ikx ■ , Tx 

■/-vr/e ZTT 



(3.6) 



where u)k = (2(1 - cos{ke))/e^y/^. It can be verified that (|3.6| ) reduces to (|2.6|) in the 
limit e ^ 0. As is clear from Fig ^, the agreement between the analytical and numerical 
computations is quite satisfactory. 

Finally, we turn to the numerical results at large t. These are shown in Fig ^ A best fit 
to the data with a power law t " gave an optimum value of a = 0.61. However an equally 
good fit to the data was obtained by the function ai/ Vf + a2/t, with the second subleading 
term contributing only about a 10% correction at the largest t (we found 02/01 = 0.65). 
This second fit is consistent with diffusion — assuming this is the correct form, we obtain the 
estimate for the numerical prefactor in ( p.5|) : 

(3.7) 



B ^ 3.32. 



In their studies of the classical lattice model, Reiter and Sjolander [|^, also computed 
the spin diffusivity. Diffusion is a property of the -C 1 regime, and we do not expect 
their perturbative techniques to be exactly valid. Combining the classical values of Xu{T), 
^(T), in (|1.24| , |1-25| ) with (|33| ), their result translates into the value B = 1. This value is 
clearly inconsistent with our numerical result above. 
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IV. INTEGRABILITY AND DIFFUSION 



This section will examine an integrable toy model of spin transport; readers interested 
mainly in the experimental implications of our results so far can move ahead to Section 0. 
Others, not interested in the details of the toy model, may want to jump to Section |1V Q 
where we will discuss general implications of the toy model solution on spin transport in 
integrable systems. 

The toy model we shall introduce is a variant of an effective model, considered in Ref , 
for the dynamics in the regime T < A for gapped chains. Here, our strategy will be to 
introduce the model as worthy of study in its own right, as it is simple enough to allow 
determination of the spin density correlator in closed form at all times. In the long time 
limit, the correlator has a diffusive form, and so this example proves that there is no general 
incompatibility between integrability and diffusion. Further motivations in examining this 
model are: 

(z) We shall show that the short time behavior of the toy model is very similar to our result 
( p.6|) for the continuum wave model (|1.14|J1.18| ). This is suggestive, and indicates that the 



long time diffusive behavior in ( p.4| ) is not an unreasonable postulate, 
(u) The correlators of the integrable toy model can also be studied for the case of finite 
system of size L with periodic boundary conditions. This allows us to carefully examine the 
interplay of the limits t — oo and L ^ oo. The diffusive form only appears if the L — oo 
is taken first. In the opposite order of limits we find recurrent behavior with a great deal 
of structure dependent upon the microscopic details of the model. We think this issue of 
the orders of limit of t — > oo and L — >■ cxo is of considerable relevance to recent studies of 
T > transport in integrable systems and this will be discussed further in 

Section |VC|. 

We begin by describing the toy model. Place N point particles of equal mass at positions 
Xi (i = 0, 1, . . . A^ — 1) which are chosen independently from a uniform distribution on a circle 
of length L. Now independently give each particle a velocity Vi, drawn from some distribution 
g{v), and a 'spin' mj, drawn from some distribution h{m). As the system evolves, the 
particles will move in straight trajectories, transporting their spin along with them. This 
will happen, until two particles collide, and we now have to describe the nature of such 
collisions. We will restrict the collisions to satisfy the important constraints of conservation 
of total energy, momentum, and spin in each collision. The first two are already sufficient to 
determine the fate of the velocities (see Fig if we consistently label the particles from left 
to right [i.e. as we move around the circle anti-clockwise, we always encounter the particles 
in the order xq, Xi, X2, ■ ■ .xjy^i), then the particles will simply exchange velocities in each 
collision. In other words, in a collision between particle i and particle i + 1, the velocity 
of the particle i after the collision is that of particle i + 1 before the collision, and vice 
versa. How about the fate of the spins rrii and mj+i ? In principle, we can choose numerous 
possibilities interpolating between zero to total reflection, consistent with conservation of 
total spin J2i ^i- The exactly solvable models are the two extremes: zero or total reflection. 
The case of zero reflection is rather trivial and leads only to simple ballistic transport of 
spin along straight lines. We will therefore consider only the case of total reflection here: 
in the convention we are following of labeling the particles here, this corresponds to the 
statement that each rrii is a constant of the motion (see Fig^). So to summarize: in each 
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FIG. 8. Collision of two particles. They initially have velocities v' and spins m, m' . After 
the collision they exchange velocities but retain their spin. 

collision the particles exchange velocities but their spins 'bounce off' each other. We note 
that, although no exact solution is known, we expect the long-time correlations of a model 
with only partial reflection to be quite similar to that of the total reflection case, but with 
renormalized transport coefficients. 

We shall be interested here in computing the correlators of the 'spin' density, L{x,t), 
defined by 

Af-l 

L{x, t) = ^i^i^ - W)' (4-1) 

i=0 

where Xi(t) are the positions of the particles: these consist of piecewise straight lines which 
reflect at each collision. We shall compute the two-point correlator of L{x,t), averaged over 
the ensemble of initial conditions defined above. Further, we will choose our distributions 
g{v) and h{m) to be even i.e. on the average, the total net momentum and spin are zero. 
Then, because the initial momenta and spin are uncorrelated, we have 

{L{x,t)L{0,0)) = Y.{mimi,){6{x - Xi{t))6{xe{0)) 
= (m2)^(5(x-x,(t))5(x,(0)) 

i 

= {m^)pP{x,t), (4.2) 

where P{x, t) is the probability that a particle at x = at time t = is at the position x at 
time t, 

(m^) = ^m2/i(m), (4.3) 

m 

and p = N/ L is the density of particles. 

Most of our results will be on a particular, simple, velocity distribution, which is designed 
to mimic the properties of the continuum model ( |1.14|J1.1^ ): 

g(v) = ^ [6{v - c) + 6{v + c)] . (4.4) 
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So each particle is allowed to have only one of two velocities, ±c. This is similar to the fact 
that linear spin-waves in the continuum wave model also have velocities ±c. 

The remainder of this section will describe the computation of the function P{x, t) using 
the method of Jepsen [^. Although later more elegant solutions were put forward by 



Lebowitz and Percus for solving the model in the thermodynamic limit, we shall use the 
relatively cumbersome machinery of Jepsen because it allows us to consider finite systems. 

At time t = the A^ particles are at random positions on the ring. We shall put the 
origin of the coordinate system at the location of particle 0. The rest of the particles from 
0, 1, 2, ■ ■ ■ A^ — 1 are numbered such that the [i + l)-th particle is immediately to the right of 
the z-th one. As a particle moves with uniform velocity, in a "space-time" diagram we can 
represent its motion as a straight line which we'll call a trajectory. When two trajectories 
cross, there is a collision. The particles bounce off each other, and in effect the particles 
exchange trajectories. So at the beginning the zeroth particle starts on the zeroth trajectory 
and as this trajectory crosses others, the zeroth particle moves onto a different trajectory. 

Now define Ajk{t) to be one if the particle j is on the trajectory k at time t and zero 
otherwise. If we form an ensemble of systems, the the average {Ajk{t)) is the probability 
that the particle j is on the trajectory k at time t. A knowledge of Ajk{t) for all values of 
j, k, t constitute a full solution to the dynamics of the system. The solution is defined by 

-1 N-l 



N-l 



U 



^, / = 0,l,2,...,Ar-l; Y.^Y. 



u 1=0 

S[u, w] = e*"" when {n — 1)L < w < nL for each n 

Wkh = Xk-Xh + (vk - Vh)t. (4.5) 

We shall note some periodicity properties of S[u, w] and Ajk{t) here. By the above definition 

S[u, w + L\= e^"^[M, w], hence S[u, w + NL\ = S[u, w]. (4.6) 

Also, for the distribution (^^), noting the fact that \vk — fh| = 0, 2c, we arrive at 

A,k (t + ^) = ^ E e-'=- if S [u, WkH + NL^-^^] = A,k{t). (4.7) 

Let us define a time T = L/c which is the time required by a free particle to go once around 
the system. Every trajectory returns exactly to its starting point after this interval of time. 
It is now obvious that in a period of time NL/c = NT, each particle will return to its 
initial positions and velocities. Thus the Poincare recurrence time of the system, with the 
velocities chosen under ( [4.4| ), is of order A^, rather than being of the order of or larger. 



A. Diffusion in the Thermodynamic Limit 

Now we go to the thermodynamic limit and address the question of diffusion. The limit 
is defined such that A^ and L approach infinity while the density N/L = p being held finite. 
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The zeroth particle starts out at the origin and we ask what is the probabihty, P{y,t), that 
it is at position y at time t (both y and t are held finite as the limit L — > oo is taken). This 
can be written as 



P{y, t) = {5{y - xo{t))) = Y.{S{y - Xk - Vkt)Aok{t)), 

k 



(4.8) 



where the angular brackets an average over all possible initial ensembles of velocity and 
positions of particles, while keeping the zeroth particle at the origin. The average can be 
evaluated exactly, and there is a simple, closed form result: 



P{y,t) = -[6{y + ct) + 6{y-ct)]e-P'^'^ 



+ ^[9{y + ci)-9{y-ct)]e-^^\'\ 



c\t\ 



\J cH^ — y- 



loipJcH'-y' 



(4.9) 



Jo and Ji are the modified Bessel functions of order zero and one respectively. The resem- 
blance to the correlation function of the non-linear wave model given in (P^ is clear. The 



first term in ( [4.9|) is a delta function along the light cone, but its contribution decreases ex- 
ponentially with time. The second term lies within the light cone, and becomes increasingly 
important for large time. Also, if we take the short time limit of ( |4.9| ), we get 



P{y,t) = I [6{y + ct) + 5{y - ct)] (1 - pc\t\ + ...) + ^ [e{y + ct) - e{y - ct)] , (4.10) 



which is precisely of the form ( ^.6] ). However, unlike ( |2.6| ), we can now also study the long 
time limit analytically. We take this limit within the light cone with y ~ i/t, and then the 
asymptotic expansions of the modified Bessel functions yields 



P{y) 



exp 



4:Dt 



D = 

2p' 



(4.11) 
(4.12) 



which is the diffusive form assumed for the classical wave model in ( |3.4| ). As shown by Jepsen, 
this calculation can also be done for a general velocity distribution g{v), and provided the 
distribution is symmetric in v, we obtain ( [4.12| ) but with 



p Jo 



vg{y)dv. 



(4.13) 



B. Effect of Periodic Boundary Conditions in a Finite Geometry 

Here we go back to the finite system and think more about it. The fact that the Poincare 
recurrence time is only linear in iV, is because of the fact that the phase space becomes 
very restricted once we allow only two possible velocities. An interesting effect is that the 
recurrence time can be very different depending upon whether there are even or odd number 
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of particles in the system. The recurrence time is the lowest common multiple of T (the 
time required for trajectories to return to their initial position) and NL/2c = NT/2 (the 
time required for particles to come back to their initial trajectory). With an odd number 
of particles, N = 2p — 1, the recurrence time is NT = {2p — 1)T as mentioned above. But 
if we add one single extra particle to the system, the recurrence time almost becomes half 
because NL/2c = NT /2 = pT is an exact multiple of T. 

In fact, the recurrence time could be even smaller. If we take an even number N of 
particles and choose their velocities as ±c randomly, the most likely scenario is one where 
half of them have velocity +c, and the rest have velocities — c. In that case 

1 7V-1 

Ajk it + T) = -Y^ e^*^" n S [u, Wkh + r{vk - Vh)] (4.14) 

Amongst all the factors that we have in the product on RHS, half of them will not contribute 
any phase to the product because for them Vk — Vh = ^- The other half will contribute a 
phase of exactly exp(±22M) each. So the total phase contribution will be exp(±2wiV/2) = 1! 
Hence the most probable recurrence time for a random ensemble of even number of particles 
is T. It turns out that with reflecting hard wall boundary conditions, the recurrence time 
is 2T regardless of the initial conditions. 

Now let us see how the effect of the recurrence time might show up in the probability 
distribution of a diffusing particle. We shall essentially try to find the autocorrelation 
function for a single specific particle moving around in a ring for all times t <T = L/ c and 
t > T = L/ c. We shall again choose this particle to be the zeroth particle and at time t = 0, 
its position and velocity are Xq = 0, Vq respectively. The probability distribution is defined 
as before in ([4.8|) . Since Ajk{t) is periodic with a period of NT /2, all distribution functions 
will also be periodic with the same period. Moreover 

1 Af-l 

A,k{NT/2 -t) = -Y^ e-''^ n S [u, w^h + (vk - Vh){NT /2 - t)] 

u h=0 
1 7V-1 

= E n S ^kh - {Vk - Vf,)t] . (4.15) 

-'^ u h=0 

When we carry out the average over the initial conditions, the particles k and h will have 
the velocities +Vk, +Vh and —Vk, —Vh with equal probability. So the probability distribution 
function will satisfy 

P{y,NT/2-t) = P{y,t). (4.16) 

Thus we need to evaluate this function only for < t < NT /A. 

The details of the evaluation of P{y, t) are again relegated to the appendix. Let us 
represent the distribution function as 

P{y,t) = [5{y + ct) + 5{y - ct)]P^'\t) + [e{y + ct) - e{y - ct)]P^^\y,t). (4.17) 

We write down explicit series solutions for P^^\t) and P'-'^\0,t). These sums could not be 
evaluated in a closed analytic form. So we carried out the sums numerically for specific 
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FIG. 9. The probability distribution function P^^\t) plotted as a function of time 
< t < NT /A for A'^ = 50. The unit of time is T. The stronger peaks at multiples of 5 is 
due to the fact that 5 is a prime factor of 
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FIG. 10. The probability distribution function P^'^\t) plotted as a function of time for = 50 
for < t < NT /A. The unit of time is T. 
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FIG. 11. The probability distribution function P^^\t) plotted as a function of time for = 50 
(dashed hne) and = 51 (full line). The unit of time is T. Notice the disappearance of the peak 
at t = 1 for iV = 51. 



values of and L. P^^\t) and P^'^\0,t) are plotted in Figs ^ and [T^ respectively. For 
times t <^T and ^ 1 with N/ L fixed, these results reduce to the results ( [1.9D derived in 
the thermodynamic limit. For t > T, we see lots of complicated structures. Not all of them 
are well understood. But some of the prominent ones are easy to understand. For example 
at times t = nT, n = 1, 2, 3, ... we see peaks in P^^\ each of whose height turns out to be 
1 / y2TxN . The origin of this is very easy to understand. Out of the whole ensemble of initial 
conditions, a fraction of them will have exactly half of the particles with velocity +c and 
the others with velocity — c (with the assumption that A^ is even). Since these set of initial 
conditions have a recurrence time of only T, so they come back to the original distribution 
after a time T and hence contribute to this peak. Using the binomial distribution, it is 
easy to prove that the fraction of ensembles which have exactly half the particles with one 



velocity and the other half with the opposite velocity is exactly The extra factor 

of 1/2 comes from the fact that we have both left and right going initial conditions for our 
test particle. If we change the number of particles by one, A^ becomes odd. So no set of 
initial conditions will have a recurrence time T, and these peaks will disappear as shown in 
Fig O. 



C. Discussion 



Our study of an integrable toy model in this section has highlighted the extreme impor- 
tance of taking the limits of long time t — oo, and large system size L — oo with proper 
care. 

If we send L oo first, then we explicitly demonstrated the existence of spin diffusion 
in the subsequent long time limit. This is, of course, the correct thermodynamic limit, and 
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the presence of spin diffusion also makes physical sense: once the limit L oo has been 
taken, an infinite number of parameters are needed to specify the initial thermal state of 
the system — diffusion then arises when some local degree of freedom starts to sample an 
increasing number of these infinite number of random initial conditions with the passage of 
time. 

On the other hand, very different results were obtained for the limit t ^ oo at any fixed 
L. Here the integrable nature of the system was immediately evident, and we observed a 
bizarre set of recurrences dependent sensitively on details of the microscopic Hamiltonian. 
The Poincare recurrence time of our toy model was quite short, and this was clearly due to 
its integrability. No sign of spin diffusion was seen. 

A number of recent studies have examined the issue of T > spin transport in integrable 
quantum systems A model of particular interest has been the S = 1/2 

XXZ chain. At the SU{2) symmetric point (the XXX chain), the low T properties of this 
antiferromagnet are expected to be in the universality class of C in ( |1 . 1| ) at ^ = vr. So studies 
of the XXX model will explicate the nature of spin transport at temperatures T < Tq at 
6 = TT. This is a regime for which our paper has no results (although, spin transport in the 
low T < A regime for 6 = was studied in Ref [Q). However, we have examined the higher 
temperature regime Tq < T < T^l^ in Section [lllij| , and found that our numerical results 
are not inconsistent with the presence of spin diffusion. It would then seem natural that 
diffusion may also exist for T < Tq, although, this is, of course, not a rigorous argument. 



The latest numerical evidence for the XXX chain ||33| , |35[| seems to be consistent with the 
existence of diffusion. 

Here ,we wish to issue a small caution towards the method used to study transport in 
Refs []3T| , |32| , |33| , |35[| (this caution does not apply to Ref @). These works computed a 'stiff- 
ness', which is the co-efficient of zero frequency delta function in the frequency dependent 
conductivity. By its very construction, such a quantity is defined at = in a finite system; 
so implicitly, the limit uj ^ has been taken before the L oo limit. The considerations 
of this section make it amply clear that such a procedure is potentially dangerous. 



V. IMPLICATIONS FOR EXPERIMENTS 

First, we summarize the main theoretical results of this paper. We have shown that 
there is an intermediate temperature range over which the static and dynamic properties 
of a large class of one-dimensional Heisenberg antiferromagnets are described by the deter- 
ministic, continuum model defined by (|1.14|) and (|1.18|) . For p-leg ladders of spin 5* ions. 



this temperature range rapidly becomes quite wide as Sp increases. For 2Sp even, our uni- 
versal results hold for A < T < T^l^, where A is the ground state energy gap, and T^^^x 
is estimated in (|1.23|) , while for 2Sp odd, they hold for Tq < T < Tj^l^, where Tq is an 



energy scale measuring the strength of logarithmic corrections at the lowest T. Both Tq and 
A become exponentially small as Sp increases, and so the intermediate temperature regime 
is clearly defined. The dynamical properties of such antiferromagnets are encapsulated in 



the scaling forms ( |1.21|) , which relate them to universal functions dependent only upon two 
thermodynamic parameters: the antiferromagnetic correlation length ^{T), and the uniform 
spin susceptibility Xu{T)- We obtained information on the universal functions in Sections |I| 
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and m, including exact results on the spin- wave damping, while exact results for ^ (T) and 
Xu(T) were presented in Section |. 

We now briefly review the work in the 70 's on the dynamics of classical antiferro mag- 
netic chain (Refs p|JIO[] and references therein). As we saw in Section |, there is a window 
of temperatures T^ax < T < T^l^ (^mix "was estimated in (|1.3| )), over which ( |1.24| ) and 
( |1.25| ) are valid, where our results apply to purely classical models; so there is a common 



regime of validity between our and earlier work. These earlier classical results were all ob- 
tained on studies of lattice antiferromagnets, and all used some variant of the short-time 
moment expansion to extrapolate to the long time limit by a physically motivated ansatz 
e.g. the memory function formalism; however, there is a degree of arbitrariness in any such 
ansatz. In a regime where their correlation length ^ ^ a (where a is a lattice spacing), 
and the wavevectors fca <^ 1, their results should be described the continuum model we 
have discussed here. However, we expect their short-time methods to be exact only for 

^ 1. Indeed, our paper is the flrst to provide a proper description of the scaling struc- 
ture, along with quantitative information on the scaling functions, in the non-perturbative 
regime fc^ <^ 1. Consistent with these expectations, we saw in Section ^ that our result 
( p.20| ) for the spin- wave damping for 1, and T^^x < T < Tj^l^, was in precise agreement 

with that of Reiter and Sjolander [|ry]. It should be noted, however, that our result ( |2.19D for 
the damping remains exact over a much wider window of temperatures (Tq, A < T < T^^^), 
including when Xu and have the quantum renormalized T-dependence in (|1.9| ) and ( |1.10D . 
Moreover, we believe, despite conjectures by Reiter and Sjolander to the contrary, that the 
results in Ref [|l^ are not exact for kC, <^ 1 — we saw in Section [IIIB| that our value for the 
spin diffusion constant (assuming the existence of diffusion), D, disagreed with theirs. 

Turning to experiments, single chain {p = 1) antiferromagnets with S > 1 which have 
been studied are (CD3)4NMnCl3 (TMMC) which has S = 5/2, (CioH8N2)MnCl3 

which has S = 2, and CsVCls which has 5 = 3/2. We think it would be worthwhile 
to re-examine these materials from a modern perspective, given the numerous exact results 
that are now available. 



Among static properties, neutron scattering experiments ||38| have measured the corre- 
lation length, ^(T), and these have been compared to purely classical theories in which ^(T) 
behaves like (|1.25|) . At lower temperatures, C,{T) should exhibit the logarithmic temperature 
dependence in ( |1.1CI| ), arising from quantum fluctuations. Combined with measurements of 
the uniform susceptibility, Xu(T), a rather precise test of the quantum- renormalized static 
theory should then be possible. 

Dynamic tests of the theory have focused mainly on the linewidth of the spin-wave 
excitations in the regime kC, ^ 1. The measured linewidths have been compared PD|PT 



with the prediction of the classical theory [|T0|, which yields the result (|2.20| ). This is in 
general agreement with the theory, but a quantitative discrepancy was observed I^T] for 
S = 3/2. We think it would be useful to compare the experiments with our new exact result 
( p.l9| ), while using the actual experimentally observed values of ^ and Xu- 

We think future neutron scattering experiments should also examine the interesting 
regime fc^ <^ 1. Here we have provided new, (numerically) exact results in Section |T|. 
In particular, there is some interesting physics in the structure of the frequency-dependent 
lineshapes in Figs ^ and |, and these should be subjected to experimental tests. Also, we 
can easily generate additional universal spectra at other positions in the energy-momentum 
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space, as needed. 

The results in Section III also provide quantitative predictions for NMR experiments 
on spin chains. The nuclear relaxation rate 1/Ti is given by local, low frequency, dynamic 
structure factor of the electronic spins. This has two contributions, one from the ferromag- 
netic component given by the correlator of L, and the other from the antiferromagnetic 
component given by the correlator of n. Let us parameterize the electronic spin Sj by 



-lYSnixi) H — L(xi 



P 



where a is the lattice spacing. If we assume that 



1 



(5.1) 



(5.2) 



where F is related to the hyperfine coupling, and — > is the nuclear Larmor frequency. 
The electron spin correlator has to be evaluated in the presence of an applied magnetic 
field H, and the electron Larmor precession can usually be neglected. However, for the case 
where there is spin diffusion, as in the assumed form ( p.4| ), then this electron precession 
must be included for the Fourier transform is divergent at low frequencies. Combining ( ^.1|) 
and ( |5.2| ) with the results of Section we obtain 

1 2 r ^/mN 1/2 



1 



— = F 



AS' 



In 



A- 



MS /J L 



e(T)Xn±(T) 



$Ko) + 



TxuiT)ia/py 
V2DH 



(5.3) 



where D is estimated in (p.5|J3.7| ). If we ignore logarithmic factors, the first antiferro- 
magnetic term in ( ^.3|) is of order F/T while the second ferromagnetic term is of order 
(F/T)(Ta/c)^(T/if)^/^; either term could be dominant, depending upon the magnitude of 
H. Further, our result (|573|) has assumed the existence of spin diffusion, but we expect 
that ( ^.3|) will provide a reasonable quantitative estimate of the H and T dependence for 
experimental purposes, even if this assumption is not entirely correct in its details: there is 
clearly a long-time tail in Fig |^, even if it is not precisely diffusive. 

Finally, we compare our theoretical predictions with quantum Monte Carlo simulations 
on odd leg ladders. Numerical results for Xu{T) have been obtained recently by Frischmuth 
et al. |^2| on antiferromagnets with S" = 1/2 and p = 3,5. We compared their results for 
p = 5, with our new result ( |l.y| , |1.12| ): this is shown in Fig |r^. The fitting parameters in this 
comparison are the values of c and Tq. There is an arbitrariness in choosing the ranges of 
T over which to fit the intermediate T prediction ( |1.9D , and this can lead to some variation 
in the values of c and Tq. A reasonable set of values are c = 2.06Ja and Tq = 0.0058J, are 
used in Fig |12|. The value of c is roughly consistent with that estimated earlier |^ , but the 
value of Tq does appear to be rather small. 
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FIG. 12. Comparison of the numerical results of Ref. (plus marks) for the uniform suscep- 
tibility, Xui of a 5-leg ladder with = 1/2 with the theoretical predictions of (l.£, 1.12| ) (dashed 
line). All exchange constants are nearest neighbor and have magnitude J, and the lattice spacing 
is a. The value of Xu is per rung. 



APPENDIX A: CROSSOVER ENERGY SCALES FOR 6* = vr 



This appendix will derive the relationship ( |1.12| ) between the two energy scales To and 
^Ms associated with the quantum 0(3) non-linear sigma model at 6' = vr. This model has a 
flow from the high energy fixed point at = to the low-energy fixed point at g = Qc- The 
fiows near both fixed points are marginal: is the energy scale characterizing the fiow 

away from the the g = fixed point, while Tq is the scale characterizing the fiow into the 
g = Qc fixed point. These scales appear in the logarithmic corrections that appear in both 
the high temperature (Eqns (|1.9|Jl.lCl| )) and low temperature (Eqn ( p..7|) ) limit. A complete 
Bethe ansatz analysis of the fiow between the two fixed points matches the two scales, and 
leads to the relationship (|1.12|) . 

We will perform the matching by considering various limiting regimes of the free energy 
density, J-'q, of Zq as a function of T and H. This is will allow us to make an intricate series 
of mappings between numerous results which have appeared recently in the literature. 

First, let us consider the low T and low H regime, where T,H <^ Tq. In this regime, 
the model is in the vicinity of the g = gc fixed point, which is the k = 1, SU{2) Wess- 
Zumino-Witten model. We know from ( |1.7| ), that for i7 ^ T ^ Tq, the free energy has a 
contribution 



if2 



1 + 



1 



21n(To/r) 



+ . . . 



< T < To 



(Al) 



From this we anticipate that for T ^ i7 ^ Tq, we will have a corresponding contribution 



1 + 



2\n{C,To/H) 



+ 



(A2) 
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where Ci is a universal number we would like to determine. The universality of Ci implies 
that we can use any model which is in the vicinity of the k = 1, SU{2) Wess-Zumino- 



Witten fixed point. In particular, Lukyanov ||2^ has recently computed the detailed H 
and T dependence of the free energy of the S = 1/2 antiferromagnetic chain with nearest 
neighbor exchange which fulfills this requirement, and we can use his results to obtain Ci. 



In particular, from (3.18) and (3.20) of Ref we determine that 



Tn 



TT 



1/2 



,7+l/4j 



(A3) 



for the nearest neighbor S = 1/2 antiferromagnet, and that 



Atic 



1 + 



2 \n{27re-"'To/H) 



+ 



H<^To, T = 0. 



(A4) 



We emphasize that the result ( |A3| ) is non-universal, while ( |A4| ) is universal i.e. only the 
latter is a property of the continuum 0(3) non-linear sigma model at = tt. In the result 
( |A^ ) we have set T = 0, as this is the limit in which we shall use in the following. 

Let us now consider the vicinity of the g = fixed point for non-zero H at zero temper- 
ature i.e. for H 3> ^jjs ~ at T = 0. In this case, a renormalized perturbation theory in 
g can be used to determine the free energy density, and this was carried out by Hasenfratz 
et al [^. They obtained 
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H » Ai75> 



T = 



(A5) 



Finally, we need to match ( |A4|) and (|A5|) by computing the T = free energy for if 7^ in 
both the low and high field limits of the 0(3) non-linear sigma model at ^ = vr. Fortunately, 
precisely this computation was carried out by Fateev, Onofri and Zamolodchikov [^]. They 
obtained (|A5|) in the limit H ^ A^g^, while the result for the opposite limit H <C 
(4.97) of Ref. [|3[: 
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1 + 
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; iy<A 



MS' 



T = 0. 



(A6) 



Comparing ([A6| ) with ( |A4|) , we immediately obtain ( |1.12| ) 



APPENDIX B: SHORT TIME EXPANSION: DETAILS 



We start with the expression ( |2.2| ) and insert ( p.4| ) into it and the expand to first order 
to find 



Zc = J V'iljVip*V(f)V(f)* exp{- J dxC), 



(Bl) 



where 
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2< 



(B2) 



Here we have dropped some additive constants and integrated by parts in places. As adver- 
tised before, the magnetic field in the z direction parametrised by m? adds a mass term to 
the action and makes the (00*) propagator infrared finite. The last term in the action arises 
from the Jacobian of the delta functionals 5{n ■ L)(5(n^ — 1) in the measure. This Jacobian 
is infinite for a continuum system. To regularize it, we can introduce a discrete lattice in 
space. In that case the parameter a is the lattice constant, or equivalently the volume of 
the Brillouin zone. 



^ fcgBZ 



(B3) 



This can be finite only in a finite volume system, but we shall carry it through and ultimately 
it will cancel all ultraviolet divergences arising from unrestricted momentum sums over the 
{^jiIj*) propagator. In terms of ip and 0, the (L(a;,t) ■ L(0,0)) is written as 

(L(x, t) ■ L(0, 0)) = (^(x, t)^*(0, 0)) + c.c. 

+ -M{x, t)ij*{x, t)0(O, 0)^*(0, 0)) + c.c. 



+ ^{i:{x, t)0*(x, t)ij*{0, 0)0(0, 0)) + c.c. 



(B4) 



Using ( |2.4D in the dynamical equations (|2.1|), we find the following equations of motion: 
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dx 



(B5) 



These equations differ from ( |2.5| ) in that we have absorbed a factor of c = ^^^^ into the 
definition of time so that time and distance have the same units. 

To solve for ip{x,t), we make a Fourier transform in space of ( P5| ) and convert it to an 
initial value problem given by 



d_ 

dt 



ij{k,t) 
<j){k,t) 



-ik'^(f){k,t) 
~itp{k, t) 



+ 



A{k,t) 
B\k,t) 



(B6) 



We have here written the nonlinearities as the inhomogeneous part of a set of first order 
equations. The solution to the initial value problem is given by 
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(B7) 



The columns of the 2x2 matrix }C{t) is made up of the two hnearly independent solution 
vectors of the homogeneous problem. 
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Armed with this, we can solve for the fields in terms of the initial conditions and iterate 
the solution to go to higher orders by plugging the solution back in A and B. Also calcula- 
tions simplify very much if we work out everything back in real space. So at last we write 
down the final iterative form from which the entire perturbation series can be generated. 

t) = ^ [i){x + t, 0) + - 0)] + + t, 0) - 0(x - t, 0)] 

1 /•* 



+ - ^ dr y" dx' [5{x - x' + t - r) + 5{x - x - t + t)] A{x' , r) 

+ 1-1 dr I dx' \5{x -x' + t-T)-5{x-x' -t + r)] ^B(x', r) (B9) 
2 Jo J ox' 



(j){x,t) = -[(j){x + t,0) + (f){x-t, 
1 /■* 



dx' [e{x -x' + t)- e{x -x' - t)\ ^{x , 0) 



+ - / dT I dx' \b{x - x' + t - t) + ^{x - x' - t + t)\ B{x', t) 
2 Jo J 

-- f dr I dx \d{x ~x' + t-T)-d{x~x' -t + r)] A{x' , r) 
2 Jo J 



(BIO) 



To simplify notation let us denote the zeroth order solutions for the fields as ■?/'*^*^-*(a;, t) and 
0(0) (x,t). 



t) = ^ [^(x + t, 0) + V^(x - 0)] + + t, 0) - 0(x - t, 0)] 



1 



i d 



0(0) (x, t) = ^ [0(x + t, 0) + 0(x - t, 



2dx 



dx' [e{x -x' + t)- e{x -x' - t)] ^{x', 0) 



(Bll) 



Now we can proceed to evaluate each of the correlation functions in (|B4| ) . As an example 
let us look at (?/'(x, t)^/'*(0, 0)). The correlations to zeroth order are 



(7/'(o)(x,t)^(o)*(0,0)) 
(0(o)(x,t)0(o)*(O,O)) 
(0(o)(x,t)^(o)*(O,O)) 
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(B12) 
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To compute the one loop correlation {ilj{x,t)ilj* (0,0)), we write down 



(^(x,t)V^*(0,0)) = (V^(°)(x,t)V^W*(0,0)) 

+ ^ ^ dr y" dx' [6{x - x' + t - t) + 6{x - x' - t + t)] x 



+ ^ [ dr f dx' \6(x - x' + t - r) - 6{x - x' - t + r)] x 
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On the second and third terms of the RHS, we can replace by (jy^^^ijj^^^) since it is 

already first order in 1/,^. Now all the correlation functions on the RHS are expressed in 
terms of correlation functions at time t = 0. So we evaluate them using the partition function 
([Bll). Other correlations can be evaluated in the same manner. So we shall only write down 
the final answers here. 



(^(x, t)^*(0, 0)) = ^ [6{x + t) + 6{x - t)] (^1 - ^ j 

t)r{x, mo, o)r{o, o)) = ^ [e{x + 1) - eix - 1)\ 

(^(x, tW{x, t)r{0, 0)0(0, 0)) = [6{x + t)+ 6{x - t)] (l + . 



(B14) 
(B15) 
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Adding them all and taking the limit m — leads to the expression quoted in ( p. 61) . 
Note that all terms which diverge as m — 0, cancel each other only when we evaluate the 
0(3) invariant correlation (L(x, t) ■ L(0, 0)). 

APPENDIX C: EVALUATION OF THE ONE-PARTICLE DISTRIBUTION 

FUNCTION 



Let us start with equation ^78| . Following Jepsen's notation we write 
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with the definitions that 
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The second term in the RHS of ( |d] ) clearly represents the probability that the zeroth 
particle stays in the zeroth trajectory at time t, while the first term represents the probability 
that it has been scattered to another trajectory. The second term can be easily simplified 
using the definitions given above and in (|4.5|) . It is 



(5(y + ct) + Siy-ct))P^'\t) 



- (1 + e*P" 



2ct - pL . ■ 
1 + 

Lj 



1) 



n N-\ 



1 + e 



-tpu 



2ct - pL , 
1 + ^-^{e 



1) 



-, N-\ 



Here p is defined so that 



0<2ct~pL< L. 



(C3) 



(C4) 



can't be evaluated in a closed analytic form. But we have evaluated it numerically 
for fixed values of N and L. P^^\t) has been plotted in Fig ^ for = 50, L = 1, c = ±1. 
As mentioned before, the peaks at times T, 2T, 3T, . . . are due to the exact recurrence of 
some configurations. We plot this function again for = 50 and 51, L = l,c = ±1 in Fig 



Tl| . Note that the peaks at times T, 2T . . . disappear. 

The first term will be handled in a manner similar to above. For simplicity we shall 
evaluate only the autocorrelation part. To do this, we set y = riL where n is an integer. 



p(2)(t) = Y^{5{xk + Vkt-nL)Aok{t)). 



(C5) 



Using previous definitions we find that 
^jv — I. 

f"^ dx 



— E — \0{nL - ct) - e{nL - ct - L) + d{nL + ct) - e{nL + ct - L)] x 

4A'' „ L 



{S[u, nL + ct] + S[u, nL — ct]) 



{S[u, nL — X — ct] + S[u, nL — x + ct]) 



/o 2L 

Here it will be convenient to define an integer p as before such that 

< ct - pL < L. 
Carrying out the integrals above we get 



(C6) 



(C7) 
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X 



L-{ct- pL) 



+ 



2N~2 

ct — pL 



,i(2p+l)n 



+ C.C. 



(C8) 



The above function P^^^ (t) has been plotted in Fig |10 
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